The quadrupole moment of slowly rotating fluid balls 
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In this paper we use tlie second order formalism of Hartle to study slowly and rigidly rotating stars 
with focus on the quadrupole moment of the object. The second order field equations for the interior 
fluid are solved numerically for different classes of possible equations of state and these solutions 
are then matched to a vacuum solution that includes the general asymptotically flat axisymmetric 
metric to second order, using the Darmois-Israel procedure. For these solutions we find that the 
quadrupole moment differs from that of the Kerr metric, as has also been found for some equations of 
state in other studies. Further we consider the post-Minkowskian limit analytically. In the paper we 
also illustrate how the relativistic multipole moments can be calculated from a complex gravitational 
potential. 
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I. INTRODUCTION 

Sources for the Kerr metric have been sought for long, 
but with little success. The only known exact solution 
that can be matched to the Kerr metric seems to be the 
rotating dust disk, with maximal allowed angular mo- 
mentum, found by Neugebauer and Meinel jl| • Some nu- 
merical studies with different equations of state for the 
source find quadrupole moments differing from that of 
the Kerr metric, see e.g. 0, d, 0, [1]. In some recent 
papers [1, 0| it is shown that in the rigidly rotating case 
incompressible fluids and fluids with polytropic equation 
of state cannot be sources of the Kerr metric in the post- 
Minkowskian limit. These results are in accordance with 
the general expectation that the ellipsoidal shape of the 
rotating fluid ball produces an extra contribution to the 
quadrupole moment which should also be present in the 
corresponding quadrupole moment of the external field 
0, SL For a general review of relativistic rotating stars 
see Q. 

In the present work we use the second order formalism 
for slowly and rigidly rotating stars, developed by Hartle 
p^ . to study the quadrupole moment and its deviation 
from that of the Kerr metric for different classes of pos- 
sible equations of state. For some earlier applications of 
the formalism see, e.g., 0, H, [ll| and for a comparison 
with numerical solutions of the full Einstein equations 
see (3]. The relativistic multipole moments of the vac- 
utim exterior metric up to order two are calculated using 
an algorithm developed in [T3|. 

The field equations for the fluid region to second or- 
der in the small rotational parameter will be solved 
numerically using fourth order Runge-Kutta. When im- 
posing an equation of state this system can be rewritten 
as a first order system of ordinary differential equations 
for nine functions, but a subsystem for six of these func- 



tions will be sufficient to study for our purposes. As- 
suming regularity at the centre the solutions will depend 
on three constants of integration, corresponding to ze- 
roth order central density or pressure, the magnitude of 
the angular velocity and one more second order small 
constant. If the solutions are required to be asymptoti- 
cally flat this second order constant will be determined in 
terms of the other constants. Due to scaling invariance 
in the angular velocity it may be given a flxed value in 
the numerical runs. Hence, given an equation of state, we 
need only vary the central pressure or density when scan- 
ning the solution space. The solutions are then matched 
to a second order axisymmetric vacuum solution using 
the Darmois-Israel procedure jisl . [Tsj . This metric 
includes the general second order asymptotically flat sta- 
tionary axisymmetric vacuum solution as a special case. 

We also consider the post-Minkowskian limit analyti- 
cally by expanding the field equations in the small pa- 
rameter A = GM/ric^ and make a comparison with the 
results of 

The paper is organized as follows: In section |TT] the 
method is briefly described and the held equations are 
presented. Also the second order vacuum metric is given 
and its relativistic multipole moments up to order two 
are calculated. The matching procedure is described in 
section IIIII and the integration constants for the vacuum 
solution are solved for in terms of the values of the inte- 
rior solution on the matching surface. In section [TVl the 
equations are rewritten in a form suitable for numerical 
integration. The results of the numerical runs are given 
in section |V] and finally a post-Minkowskian analysis is 
made in section IVTl 
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To second order the metric of a slowly rotating ax- 
isymmetric object, both in the interior fiuid region and 
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the outside vacuum region, can be written as 

ds^ = {l + 2h)A^dt^ - {l + 2m)^dr^ - 

(1 + 2k)r^ \de^ + sin^ 9 {dip - ujdtf] , (1) 

where uj is first order and h, m are k are second order 
in the rotational parameter [lo| . The requirements of 
regularity at the centre and asymptotic flatness imply 
that the first order function uj depends on r only. The 
second order functions h, m, k can be given as 



h = hn + 



m 
k 



+ h2P2{cosO) 

Too + TO2P2(C0S ( 



= fc2P2(c0S( 



(2) 



where /iqjTOo and /i2,TO2,A:2 are functions of r only and 
P2{x) = ^{ix^ — 1) is the second order Legendre poly- 
nomial. This result follows from reflection symmetry in 
the equatorial plane, from that the equations for /i, to 
and k separate with the ansatze h — X^i^o ^i('')^j(cos6') 
etc., and from the fact that there are no inhomogeneous 
terms containing lo in the equations for hi , ki and rrii for 
i> 2. For more details see jlOf . 

The matching of the two spacetime regions happens 
via the application of a coordinate transformation 
if + U,t 'm the fluid region. It means that the inner fluid 
region rotates with respect to the distant stationary ob- 
servers with angular velocity O. This parameter 17 is 
considered to be the small expansion parameter with re- 
spect to which (jj is first order and the other corrections 
h, TO, k are second order. In addition to this, we can also 
rescale the interior time coordinate first by a zeroth or- 
der constant C4, and then later by a second order small 
constant C3, i.e. t ^ Ci{\ + c^)t. 

The matter content of the interior is modelled by a 
perfect fluid 



Tab = {P+ p)UaUb - PQab 



(3) 



The coordinate system used in ([T]) is assumed to be co- 
moving with the fluid, i.e. the 4-velocity is assumed to 
possess the form 

u'^ = 0, 0, 0) = ((1 - h)/A, 0, 0, 0) (4) 

which also implies that the shear of the fluid is zero so it 
rotates rigidly. 



A. The field equations 

In this subsection we list the field equations relevant 
to various orders. A similar system of equations with a 
slightly different choice of variables was given in [l^. If 
no equation of state is specified then the only equation 
one gets to zeroth order of the rotational parameter is 



the pressure isotropy condition 1 
as 



2 which reads 



B 



d^A dirA) d{B/r) 



dr"^ 



dr 



dr 



A 

7^ 







(5) 



Making use of G° = P% and GS = the energy 
density and pressure of the non-rotating configuration 
reads as 



1 



d{rE 



2M 



dr 

B^dijA^ 
dr 



1 



(6) 



To first order in the rotation parameter the only relation 
follows from G'^ = 



dr 



^ r^^^ +^r'uj— - =0 



A dr J 



dr \ A 



B 



(7) 



The second order Einstein equations yield the following 
conditions. From G^ 2 = one gets 

r-^ (ft.2 + k2) +r {h2 - TO2) - ^2 - TO2 = . (8) 

The pressure isotropy condition in the angular directions, 
G'^2^G^3, gives 

aB'^ fdojV „ ^B d f B\ „ 

(9) 

The equality of the pressure in the angular and radial di- 
rections, i.e. G^ 1 = G^ 2 gives two equations. After ehm- 
inating the derivative of h2 using ([8|) one obtains from the 
P2(cos6') part 

B^dA f dk2 \ 7 , d [B 

m2-Ak2-^h2-\r^^(^^ =0, (10) 

while the ^-independent part gives an equation for toq 
and ft-o- 

The energy density function can be decomposed as p = 
Pn + P2, where p2 = P20 + P22P2{cos9) and p20 and P22 
are second order small functions of the coordinate r. The 
analogous decomposition of the pressure is similarly given 
by p = po +P2 = Po +P20 +P22P2 (cos 9) . We get one more 
equation for toq and ho if we assume that the equation of 
statep — p{p) is unchanged to second order. For details 
see [nl. 



B. Equations of state 

To complete the system one more equation, e.g. an 
equation of state 



P = pip) , 



(11) 
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has to be specified. In this paper we will consider equa- 
tions of state in the form: 



p — dip + d2 



Pc 



1/7 



(12) 



Here pc is the central pressure. Equation (|T2|l includes 
some of the more used approximations for the equation 
of state of dense stars, like Newtonian polytropes, rela- 
tivistic polytropes with di — 1/(7 — 1), linear ones as 
well as the incompressible case. 

We use units such that the speed of light c = 1 and 
SttG = 1. If then the unit of length is taken as a meter, 
the relation between our units and SI units are given by 



Psi 



5.358 X lO^S-^kg/m^ 



xsi = 4.813 X 10'*^— kgm^^s" 
disi = di 



(13) 



where x = pc^ ,p, d2 or ^3. 

Note that mo and /lo do not appear in equations l[5|). 
0, ®, dl]), (dni) and (HH) and hence this subsystem for 
A, B, uj, 1712, ^2 and h2 decouples. Notice also that the 
equations can be solved order by order and that the equa- 
tions contain m2 only algebraically. In section |V] the sys- 
tem will be reformulated as a coupled system of six first 
order ordinary differential equations. Due to the require- 
ment of a regular centre the solutions to this subsystem 
will only depend on three constants of integration. 



ho 
h2 



1 



= -mo = 



C2 



r - 2M 
3cir {2M - r) log 

+2ci — (3r^ - 6Mr - 2M^) 

r ^ ' 2M -r 

2M\ , 

— j-v 



2M 
r J r 
r- M 



+ {M + r) 



2M 



m2 



3ci{r^ - 2Af2)log 1 



-2ci — {2M^ - 3Mr - 3r^) 
r 

6a — j /i2 . 



a2^(2M + r) 



[2M^ - r^) qi 



(14) 



In this approximation, the slowly rotating solution is 
characterized by the mass M, the first order small rota- 
tion parameter a, and the second order small constants 
ci, C2 and qi. When qi takes the value zero the metric 
is known to be the general asymptotically flat stationary 
and axisymmetric vacuum metric to second order (see 
e.g. [3). It can be easily checked that the solution is 
of Petrov type D only if both ci and qi are zero. The 
metric is then equivalent to the Kerr metric to second 
order with mass M ^ M — C2- 

When qi ^ the metric cannot be asymptotically flat. 
It is important to keep in mind, however, that without 
the inclusion of this constant the matching conditions on 
the zero pressure surface are overdetermined in general 



C. Vacuum metric 

In the exterior vacuum region we will use a frame 
adapted to the asymptotically non-rotating observer. 
Solving the field equations detailed in Section III Al by 
imposing p = p = the metric functions for the vacuum 
region are given as follows [13, [l3| ^ 

=B^ = 1- 2M/r 



_ 2aM 



To make the norm of the timelike KilUn g fie ld equal to unity 
at spatial infinity, the expression of ho in [131 has been slightly 
modified through a second order coordinate transformation of 
the time coordinate. 



D. Multipole moments of the vacuum metric 

The notion of relativistic gravitational multipole mo- 
ments for static asymptotically flat vacuum spacetimes 
was developed in I20| bv Geroch and later extended to the 
stationary case in [2ll| and [13] by Hansen and Thorne. 
These moments are defined on the 3-space of the timelike 
Killing trajectories. 

A 3-space {A4, h) with positive definite metric h is said 
to be asymptotically flat if it can be conformally mapped 
to a manifold ( , h) with the following properties 

(i) A^ = A4 U A, where A is a single point, 

(ii) - = 0, b,b,n\^ = h,\^ 

where /ly = Q?hij. 

From the timelike Killing vector field K'^ one con- 
structs the two scalar functions / = K'^Ka and t/i, 
where the later is obtained from the curl of K°' , V,a — 
^abcd.K'^ K'^ '^ , that is a gradient due to the vacuum equa- 
tions. The complex gravitational potential then reads 
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in terms of the Ernst potential £ = f + itp. It is given the 
conformal weight —1/2, so that f = fl^^^^^. For axisym- 
metric spacetimes the metric is completely determined 
by^the value of the potential on the axis of symmetry 

The multipole tensors on the 3-space of timelike Killing 
trajectories, with coordinates a;*, i = 1,2,3 and metric 
given by the projection operator hab = -fgab + KaKt, 
are then defined recursively as 

(1)^ 



(n) 



.fe„> 



-n(2n-l)i?<fe,fe,P^ 



(n-l) 

k„ + i> 



(16) 



where < ki...kn+i > denotes the symmetric and trace- 
free part. Di and Rij are the covariant derivative and 
Ricci tensor, respectively, with respect to the 3-metric hij 

p. The are defined correspondingly in terms of 

the tilded quantities. 

For axisymmetric spacetimes the multipole moments 
are given entirely in terms of the scalar moments defined 
as 



1 ~ 



in) 



1^ ki...k„' 



(17) 



in terms of the axis vector n* [2l| . 

In, e.g., [l2| and [l^l algorithms for calculating the 
multipole moments of a stationary axisymmetric space- 
time are developed. We here use the method of [l^ to 
calculate the moments up to second order. First the met- 
ric is transformed to the canonical form 

^ _ ^^^)2 _ j-1 ^g27 (-^^2 _^ ^^2^ ^ ^2^^2j ^ 

(18) 

where /, 7 and iu are functions of p and z, through the 
coordinate transformation 



p = Ar sin 6 1 
z = (r-M) 
+r'A^ ( ho 



h-i + k2 - - sin^ 



1 + 2ho - - sin^ ( 



ih2 



9{h2 + k2) 

2k2 - m2 



■ Sin^ 6 {h2,r + k2,r) 



cose (19) 



that holds to second order. 

From the timelike Killing vector if" = 5^ we find the 
two potentials 



/ = K'^Ka - .goo - (1 + 2h)A^ 



2 2-2 

r Lo sm ( 



and 



2aM 



■ COS0, 



(20) 



(21) 



and hence ^ = (1/f)^. As shown in [l2|, the first (up to 
n = 3) scalar moments P„ are given by the coefficients 
771,1 in the expansion 



I(P = 0) 



■^n=0"''n^ 



(22) 



of ^ on the axis, that in the original coordinates corre- 
sponds to = (or 7r). Expansion of ^ on the axis in 
the original coordinates gives to second order in the ro- 
tational parameter 



M 



Ma 



2„2 



r~M (r - M)2 r{r - Mf 
2 

(/7o + h2)A^ ' 



r-M 



(23) 



From one finds that r = z + M to zeroth order, 
which is sufficient for the three last terms in ((23)) since 
they are already first and second order. The first term 
need to be expanded to second order in the rotational 
parameter, giving 



M 



M 



r-M 



z - 2ho{r - M) ~ r^A'^ho.r 

— + ^[2ho{r-M)+r^A^ho,r] (24) 



where once again r = z + M can be used in all second 
order functions. Finally, expressing ^ in z, one obtains 
to second order 

I = M - C2 iMaz - M (^a^ + ^M^a^ z^. (25) 

Hence the mass is given by Af — C2, the angular mo- 
mentum by J = Ma and the quadrupole moment by 
Q = -M (a^ + ^M^ci) . We will be interested in the 
relative deviation of the quadrupole moment from that 
of the Kerr metric 



AQ _Q- Qkc 



leM^ci 



(26) 



An expansion of the exterior metric for large r gives 
the following leading terms of goo (with qi — 0) 



500 



^ 2M (1 - U) _^ 2MP2{cos0) {a^ + fM^a 



(27) 

i.e., the associated Newtonian quadrupole moment reads 
as (cf., e.g. M) 

Qii = Q22 - -Q33/2 = 2M (a^ + yM^ci") (28) 



obtained from the curl of if. 

Introduce the coordinates p = p/{p'^ + z'^), z = z/{p'^' 



z ). A suitable conformal factor is then 



P 



in an asymptotically Cartesian system with the 3-axis 
along the axis of rotation. Hence it is, up to a factor of 
4, the same as the relativistic moment. 
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are 
and 



III. MATCHING 

We here briefly describe the matching procedure, in 
which the Darmois-Israel junction conditions 
used. For more details the reader is referred to 
also to [iBl for a general discussion on the matching of 
axisymmetric bodies in second order perturbation theory. 
In the fluid region, the matching surface S is defined by 
the condition of vanishing pressure, p — 0. In the limit 
of no rotation, the matching surface is the 3-dimcnsional 
cylinder r — ri. For slow rotation the equation of the 
matching surface S is 



(29) 



with 



P20 +P22^2(C0SI?) 



dpo 



= ^0+6^2(0081^) 



I r=ri 



where the constants ^ and ^ are given by 
1 



(30) 



Co 



12r5^f (4) 



l2rA^^ - 12mo^ (rA^) + rM ^ 

riry /-try " ' \ /-ft-< 



dr 



dr 



duj \ ' 



dr J 



and 



3A 



dA 
dr 



(31) 



In the vacuum exterior region suitable hypersurfaces 
for matching are determined by 



r = ri + X = xo + X2-P2(cos?9), 



(32) 



where Xq ^-nd X2 are constants to be determined by the 
matching conditions [265 • 

In order to find isometric embeddings of the matching 
surface S in the vacuum and fluid domains, we equate 
with each other the respective induced metrics , ds^ 

and induced extrinsic curvatures K^^\ K 



K(^^s^K\s (33) 

where K is defined by 

K = Kabdx^dx^ = /i//i/n(c;d)dx°da;'' (34) 

in terms of the unit normal Ua of the matching surface S 
and the projection operator hj' = naii^ + ^a- To adjust 
the different coordinate systems with each other we apply 
a rigid rotation in the fluid region by setting ip ^ ip + flt 
where 17 is a constant. Then we re-scale the interior time 
coordinate t € R hy t €4(1 + 0^)1 with further zeroth 
and second order constants C4 and C3 to be determined 
from the matching conditions. 



From the zeroth order matching conditions we get the 
following relations (all functions here and in the follow- 
ing first and second order relations are evaluated at the 
zeroth order matching surface r — ri): 



M=iri(l-S2) 



B 



dr 

To first order we solve for a and il 
Brl duj 



3A(B2 - 1) dr 



_ ri duo 
3 dr 



(35) 
(36) 

(37) 



From the second order equations we can solve for ci, C2, 
C3, 9i, Xo and X2 as 



B' 



Cl 



Co 



\B'^(B^ 



9r2A2(B2 _ 1)6 
-^36^2/12(1 - B^) + 12A?B'^{h2 + k2) 



(38) 



B' 



C3 = 



dB\ , r\ B^ /dcjV 

(39) 
(40) 



dLo \ ' 
^ + 



C2 



36 A2 \dr J nB^ 



- ho 



91 



^^^^^,^^^{18.2^^(5^ -1)(B^-8B^ + 1) 

-^216^^52 InS [2B'^{h2 + k2) + (1 - B^)h2] + 
36h2A^B^{B^ - 1){B^ + B^ - 8) + 

r\ B^ [(i?2 - 1)(2 -t- 11^2 - IB^)^ 

6B2(B'*-3)lnS]| (41) 

Xo = Co and X2 = 6 , (42) 

where and ^2 are obtained from pTjl and PT|) . Note 
that we did not impose the equation of state p+|) when 
calculating the matching conditions. Hence an appro- 
priate matching can be done, i.e. the vacuum metric in 
section III Cl is general enough for describing the exterior 
of any axisymmetric rigidly rotating perfect fluid ball up 
to second order. 



IV. EQUATIONS AND BOUNDARY VALUES 

In this section we provide a reformulation of the field 
equations to a form more suitable for numerical integra- 
tion. By doing this we can get higher precision at the 
origin where apparent singularities arise, moreover the 
freely specifiable constants are identified more easily this 
way. 
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A. Integrating the zeroth order field equation 

In order to simplify ([5]) it is convenient to redefine the 
functions A, h, m and k in terms of the function v as 



^ = e^ h = he 



TO = me 



k^ke'^". (43) 



The equations simphfy considerably due to the fact that 
only the derivative of v will appear. Hence we introduce 
the function z by 



z du 



(44) 



Then the zeroth order equation ([5]) becomes first order 
in z and algebraic in B [27[ 



Br 



dz 
dr 



+ ^2 _ ^Bz +1 = 



(45) 



furthermore, the pressure of the non-rotating configura- 
tion IHl) takes the form 



po = ^ (iBz -B^-1) 



(46) 



and finally to fourth order 
62 



Z2 



4 + ^« + »5-4M0| 



2zi 



bizi 



hi 



1,(2) _ 

^2 — 

(2) 

^2 = 



{izl - 16zi5i + 8bl - 105: 



42 (zi-61) 



-|^(mi-5zi) 



^(zi-35i) + ^(6i-zi)-4'\ 
2 



zi - 36i) 



3 

^ (z2 - 36zi6i + 7462 _ ^q^^^ 



The central density and pressure are given by 
Pqc = -65i , pqc = 2zi . 



(50) 



(51) 



This shows that for realistic configurations 61 < and 
zi > 0, consequently the zi — bi term in the denominator 
of /ij^^ is nonvanishing. 



B. Series expansion around a regular centre 

For sufhciently regular configurations close to the cen- 
tre the metric coefficients can be given as power series 
in r. Assuming that the central pressure and density are 
finite it follows that B{0) = z{0) = 1. The assumption of 
smoothness of the configurations at the symmetry centre, 
in the spacetime sense, implies that the odd coefficients 
in the expansions of the basic variables are zero. 

Hence, assuming a smooth centre in the spacetime 
sense, the odd powers will be omitted hereafter. Plugging 
the expressions 



B 

z 

h2 

TO2 

k 



= l + bir^ +b2r^ + ... 
= 1 + zir^ + Z2r'^ + ... 
= ujo + oJir'^ + iLJ2r'^ + ■■■ 



(1)^2 ._(2) 4 



toX -|- toX r +712 



JO) 



(1)^2 , , (2) 4 



fca 'r 



fc^ 'r\. 



(47) 



into the field equations together with the equation of 
state dTTI) justifies then that all coefficients can be given 
in terms of zi, loq and hi = h''^^ . 
To zeroth order one obtains 



then to second order 



,(0) 



k 



(0) 



0, 



(48) 



TO 



(1) 



k 



(1) 



-h 



(1) 



rWo(zi - 36i) , 



-hi , 
bi 



1 

-gP|22i 



(49) 



C. System of differential equations 

Motivated by the results of the previous section it is 
advantageous to define the new dependent variables /?, ^, 
uj, ijj, h, k and to through 

B = l + r/S , z = l + rC, 

to = LOo + roj , LU^r — U) + UJ , (52) 

h2 = rh , k2 — r{r'^k — h) , TO2 = r{r'^rh — h) . 

The closed subsystem of equations ([5|), ©, jS]), jS]), 
PH)) and pT|) then gives six first order differential equa- 
tions for the quantities /?, Q, Co, cD, h and fc, while rii can 
be solved for algebraically. The equations for the zeroth 
order quantities are given by 



dc 

dr 
dp 

dr 



1 

Vb 

1 



2rB 



{-( - 3rC/3 + rC^ + 2r/3^) (53) 
rpo + 4/3 + 3r/3^) (54) 



where the density po is given by equation of state po — 
p{po) with the pressure given by 



1 



Po = -(2C + 2rC/5-r/32) 



(55) 



The equations for the first order quantities lu and lj are 

diJ LO 

dr r 

did 1 [ ,. . . df3 



(56) 



1 

dr rB 

+ {Url3 + 4:-5rC)(:} + {7rP + 5-rC)i:j (57) 



r (4cjo + 5rc;) + rtj) -- + 4 (2/3 - () loq 
dr 
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whereas the ones for the second order quantities h and k 
are long and hence are not given here. The corresponding 
homogeneous equations, given by putting w = in ([8]), 
® and Uni), are 



dhh 2rkh — pfih hh dp 

'dr ~ r{l3-C)B ^ /3- C"*- 

dfch _ rfch (2rC - 5r/3 - 3) + 2/ih (/? - C) 
dr r^B 



(58) 



(59) 



These two last equations will be used when adjusting the 
parameters so that the solutions become asymptotically 
flat. 

Boundary conditions at r = are given as 



dp 
dr 
dh 
dr 



, dC duj du; 2 



dk Lul 



36i 



^1 (U 



zi) . (60) 



The constants bi — — poc/6 and zi — Pqc/2 are related 
through the equation of state pip , implying that we have 
three independent constants of integration (apart from 
possible constants in the equation of state). One of this 
{hi) will be determined through the matching conditions 
on the zero pressure surface in terms of the others when 
the requirement of asymptotic flatness is imposed. 

There is a rescaling associated to the rescaling of the 
rotational parameter wqi following the rule ujo — > ^(^o^ 
which induces the transformation 

(I), cD ^ 70;, 7(2;, ^Cj h,rh, k ^ 7^^, 7^^; 7^^j 
/3,C,r^/3,C,r. (61) 

Due to this scale invariance of the equations the constant 
luq can be fixed. All other configurations with a given 
equation of state and specified central pressure can be 
obtained by rescaling. 

The rescaling of the radial coordinate r — > ar induces 
a change of the parameters in the equation of state, (fT2ll . 
given by 

di^di, d2,d3,Poc ^ ^ ■ (62) 

a a a 

The dependent variables then scale as 



/3,C,w,a),m,fc 
uj, luq, ft, — > lj, ojq, ah 



P ( oj u} rh k 



(63) 



Using this scale invariance we can set one of the constants 
in the equation of state to a fixed value. 



D. Asymptotically flat solutions 

Asymptotically flat solutions can be obtained by con- 
sidering suitable linear combinations of homogeneous and 



particular solutions for /12 and k2, or equivalently for 
h and k. These are solutions of the linear equations 
([55)1 and ([55)1 and their corresponding inhomogeneous 
equations respectively. Since we solve the equations or- 
der by order, and hence k and h do not appear in the 
lower order equations, we can add any homogeneous so- 
lution to a given particular solution, i.e., k = kp + Ckh, 
h — hp + Chh, where C is an arbitrary constant and the 
subscripts p and h refer to particular and homogeneous 
solutions respectively. A solution to the field equations 
is asymptotically flat iff qi = 0. Now the expression ([^T|) 
is also linear in k2 and ^2, or equivalently in k and h, it 
has a structure like 



qi = ai + a2k + ash. 



(64) 



Hence 



qi = ai + a2ikp + Ckh) + a3{hp + Chh) ^ 
a\ -f a2kp -I- a-ihp + C{a2kh + a^hh) 
qip + Cqih 



(65) 



and qi = can be obtained by choosing C = —qip/qih- 
Also the constant ci that is used to calculate the 

quadrupole moment, has a similar structure, and hence 
the resulting ci is given by ci = cip + Ccih- 



V. NUMERICAL SOLUTIONS 

In this section we consider solutions with equations 
of state (fT2|) . with special attention to the quadrupole 
moments and their deviation from that of the Kerr met- 
ric. In particular we determine the quantity AQ/Q de- 
flned in (|26|) . The system was solved using fourth or- 
der Runge-Kutta. In scanning the parameter space, due 
to the scaling invariance, without loss of generality, we 
flxed luq = 0-1 and the method for obtaining asymptot- 
ically flat solutions, described in section PV Dl below, in- 
directly fixes the value of hi. Hence, for a given equation 
of state, we only need to vary the central pressure, 2zi. 
When considering different equations of state we use the 
scale invariance and to fix one of the constants 
d2 or o?3. When performing the corresponding rescaling 
r — > ar the quantity AQ/Q is invariant. 

The integrations were carried out until the zero pres- 
sure surface was reached. The accuracy of the code was 
checked for the Wahlquist solution [2^ . 



A. Incompressible fluids (interior Schwarzschild) 

This case was already considered in 3]. The numer- 
ically calculated quantity AQ/Q, defined in (^5]) . as a 
function of the central pressure is depicted in Figure [1] 
In Figure [2] the same quantity is given as function of the 
radius of the corresponding non-rotating configuration. 
When the central pressure approaches infinity and the 
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state p = Sp + j-iQ. Its metric is 



ds^ = fodr 



where 



/o = l 



/o (1 - «VorV2) 



{de"^ + sin^ 9d(p'^) , 



1 — arcsin ( K\ / — r 



and the pressure is 



P = (l - kVo) 



(67) 



(68) 



FIG. 1: The quantity AQ/Q as function of central pressure 
for a sequence of of fluid balls with constant density p = 1. 



radius the Buchdahl limit then AQ/Q approaches 
0.0212. Due to the scaling invariance, ([5^ and (pH)) . this 
Umiting value is independent of the density. 



o 
3 

< 




0.01 



1.2 1.4 1.6 1.8 2 2.2 2.4 
in units of 2M 



FIG. 2: The quantity AQ/Q as function of radius r\/2M for 
a sequence of of fluid balls with constant density p = \. 



B. Linear equations of state 



Linear equations of state 



p = dip + ds 



(66) 



are of interest not only because they could approximate 
some compact objects, but also since in this class one 
might have a chance of finding exact solutions, at least 
to zeroth order. 

An exact spherically symmetric perfect fluid solution 
is given by the Whittaker metric [2^, with equation of 



It is the non-rotating limit of the Wahlquist metric [28| . 
that cannot be matched to an asymptotically flat vacuum 
region [13, HO]- However, there is another rotating gen- 
eralisation of Whittaker that can be matched, although 
its closed- form is not known [3l| . 

In the limit when the central pressure goes towards 
Po/2, corresponding to k — > 0, the Whittaker metric 
approaches the anti-de Sitter spacetime and the radius 
ri blows up to infinity as y^2/ po/ k, but ri/2M goes as 
l-|-4K^/7r^, i.e. the surface approaches the event horizon. 
This is not in conflict with the Buchdahl limit f33| , since 
this spacetime does not meet the usual physical require- 
ments. Even if the these are weaker than in Buchdahl's 
original work, see e.g. [11], the central density becomes 
negative in the anti-de Sitter limit. For the correspond- 
ing second order rotating and asymptotically flat config- 
uration we find that the quadrupolc moment approaches 
that of the Kerr metric in this limit. 

Similar results hold for all configurations with di < — 1. 
For all these poc = —Poc when poc = d'i/{—di — 1), 
and since the anti-de Sitter spacetime is a solution to 
our zeroth order equations and the solution is uniquely 
given by central pressure and density this is the result- 
ing spacetime. Since the zero-pressure surface is further 
and further pushed outward when approaching the limit 
Poc = dz/{—di — 1), all these configurations also become 
infinite in extension. The ratio r\/2M once again tends 
towards one, and we also reobtain that the quadrupole 
moment approaches that for Kerr. Even if these space- 
times are quite unphysical, a study of them still is helpful 
in understanding which conditions are needed for a suc- 
cessful matching to the Kerr metric. 

When performing the numerical runs, due to the scal- 
ing invariance, ((62|) and ([63|l . we put = 1. Negative 
values of dz are excluded if the surface density should 
be larger than or equal to zero, and configurations with 
da = are not finite in size [3J|. The quantity AQ/Q 
for a sequence of values for d\ are given in Figure [31 As 
seen, for more realistic configurations with c?i > 1 it dif- 
fers significantly from zero. 
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FIG. 3: The quantity AQ/Q as a function of central pressure 
for a sequence of fluid balls with linear equation of state p — 
dip + 1. 



FIG. 4: The quantity AQ/Q as a function of central pressure 
for a sequence of Newtonian polytopes with equation of state 



, 1/7 



C. Polytropes 

1. Newtonian polytropes 

In a recent paper f?\ it was shown that slowly and 
rigidly rotating polytropes cannot be sources of the Kerr 
metric in the post-Minkowskian limit. Here we numeri- 
cally find that slowly and rigidly rotating polytropes with 
arbitrary strength of the gravitational field cannot be 
matched to Kerr either. 

/ \ 1/7 

We use the equation of state p = d2 ij- j obtained 

by putting c?i = c?3 = in (fT2|) . Due to the scaling invari- 
ance, §^ and we put da = 1. In Figure H AQ/Q is 
shown as a function of the central pressure for a sequence 
of values of 7, including the physically interesting cases 
7 — 4/3 and 7 = 5/3. As shown in [sj] the configu- 
rations are finite in size for 7 > 1.2295 (I/7 < 0.7695) 
and infinite for 7 < 1.2 (I/7 > 0.833). In the interval 
in between the central pressure determines whether the 
configuration is finite or not. As seen AQ/Q never seems 
to approach zero. Furthermore, for physically reasonable 
configurations with Pc < Pc = 1, the quantity differs sig- 
nificantly from zero. 



2. Relativistic polytropes 

For a fluid with one type of constituent particles, a 
relativistic polytrope is given hy p = CnT' , where C is 
a constant, in terms of the particle density n and the 
polytropic index 7. This equation is suitable, e.g., to de- 
scribe an ideal degenerate neutron gas. Using the energy 
conservation equation for a perfect fluid it is then easy 
to show that the equation of state becomes 



1 



in terms of the pressure p. As seen the equation of state 
approaches a linear equation of state for large pressures 
and a Newtonian polytrope for low pressures. For a dis- 
cussion of relativistic polytropes see. e.g., [IHl. The re- 
sults of the numerical runs are given in Figure [5l As 
expected the results are similar to those of the Newto- 
nian polytropes. 




FIG. 5: The quantity AQ/Q as a function of central pressure 
for a sequence of relativistic polytropes with equation of state 



1/7 



p 



1 



P + d2P 



1/7 



(69) 



D. Combined linear and polytropic 

When adding a constant to the Newtonian polytropes 
in Figure [5] typically the local minima and maxima get 
less pronounced or disappear and for large central pres- 
sures the quadrupole moments get closer to that of the 
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incompressible case, but never smaller. For combined lin- 
ear and polytropic equations of state the curves become 
more similar to those of the relativistic polytropes in Fig- 
ure O and when there is additative constant the limiting 
value of the quadrupole moment for large pressures once 
again gets closer to that of the incompressible case. 



VI. THE POST-MINKOWSKIAN LIMIT 

Since we cannot solve the equations in the fluid region 
analytically, we look at the weak gravity limit by making 
an expansion in the small parameter A — M/ri, or in SI 
units GM /ri(? . Here M is the mass of the fluid ball and 
ri is its radius. A similar study, using global harmonic co- 
ordinates and the Lichnerowicz matching conditions [3^ , 
was performed in Q and 01 ■ 

From an expansion of the equations in this parameter 
it turns out that the pressure is one order higher in A than 
the density. This can also be understood in a Newtonian 
context where the virial theorem applied to a spherically 
symmetric configuration gives {p) — kX{p), where fc is a 
constant of order unity and (p) and (p) are the average 
values of pressure and energy density. By considering the 
balance between gravitational and centrifugal forces on 
the equator of a rotating Newtonian fluid ball one obtains 
for the angular velocity that lu < \/X/ri, so it seems 
reasonable to let lo go as VA. This is also consistent with 
that LU or its derivatives only appear quadratically in Q 
and (fTO|. 

We first consider the post-Minkowskian expansion of 
the spherically symmetric system. The functions A and 
B are expanded as 



A^l + Ai + A2 + . 



B^l + Bi+B2 + ... , (70) 



where subscripts refer to the order in A. The freedom of 
rescaling the time coordinate was used to put the zeroth 
order term in A to one. We do a similar expansion of the 
pressure p = pi + p2 + ■■■ and density p = pi + p2 + ■■■■ 
Using the field equation ^ we first obtain that pi is a 
constant. We put this constant to zero so that a zero 
pressure surface can be obtained. Then 

dAi 
dr 

The first nonvanishing terms of the density and pressure 
are 



B, = -r- 



(71) 



Pi 



P2 



2 d 
dr 

Pc - 



, dAi 
dr 



(72) 



To complete the system we have to impose an equation 
of state p = p{p) , giving an equation for Ai . The differ- 
entiated version of this equation reads as 



dp{pi) dpi 
dpi dr 



-pi- 



dAi 
dr 



(74) 



where pi is given by (|72|) . We note that the total mass 
M of the configuration will be of the order A, while the 
radius ri will be of zeroth order. 

The function cj, which is first order in the rotational 
parameter 17, has the post-Minkowskian expansion 



1^1/2 + '^3/2 + 



(75) 



where the indices correspond to the appropriate frac- 
tional orders of the quantities in A. From the field equa- 
tion ([7]) we obtain that wi/2 is a constant. This corre- 
sponds to a rigid rotation of the system. It is necessary 
to keep uji/2 nonzero, since this will provide the nonlinear 
source term in the higher order equations. 

The quantities second order in the rotational parame- 
ter, i.e. /i2, m2 and k2, are expanded as 



h2 = h2l + /l22 + 
TO2 = r7l21 + m22 
k2 = k2l + k22 + 



(76) 



where the second index at each new quantity indicate 
the order in the post-Minkowskian parameter A. To first 
order in A equations ([5]), (O and PU)) yield m2i — k2i = 
— ft,2i and 



2dh 



21 



dr'^ 

where 
dpi 



dr 



dpi 
dAi 
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h2i = 



■^1/2 



dpi 



dAi' 
(77) 



_^/Ai 
dAi dr^ I 



'dAi 
dr 



A d^Ai 
r dr^ , 



'i^. (7S, 

dr 



We note that a particular solution to cq. ([77)1 is given by 
h2ip = -ujj^^r'^/3. 

Using the condition that qi, given by Equation (I4ip . 
should vanish (to first order in A) for an asymptotically 
flat solution and that the corresponding ci is given by 
ci = cip ~ cihqip/qih as described in Sec. IIVDI one ob- 
tains 



ci = - 



5^1/2 


37 log ( 






f loj 


l{rh2iH/^)U.^ 



(79) 

Here h2ih is a regular nonzero solution of the homoge- 
neous version of ([77| . Note that in general ci goes to 
infinity as l/A"*. 

From the matching condition ((35l) it follows that 



dAi 
dr 



M 



(80) 



where we assume that the mass M is positive. In [11] it 
is shown that the rotating fluid ball is oblate in shape if 



k2\r=r, H < 

n 



(81) 
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where ^2 is given by (|3ip . Substitution of fc2i = — = 
■^1/2^^/"^" ^2ih in (El) together with ([50)1 gives that the 
fluid ball is oblate iff /i2i/i|r=ri > to first order in A. 

For the incompressible case the equations can be in- 
tegrated completely and give for the asymptotically flat 
case 



h2i 



Pir 
12 



P2 ^Pc 



12 



(82) 



The integration constant ao can be transformed away by 
a first order (in A) rescaling of the time coordinate. The 
choice flo — —pirf/A makes the constant C4 in ()35|) equal 
to unity. Since /i2i/i|r=ri = 5u;^^2^i/6 > the ball is 
oblate. From this solution the radius ri of the fluid ball 
and the parameters M and a of the exterior metric are 
given by 



ri 



Pi 



M 



4(3p, 



^3/2 



3 Pi 



24 



(83) 

respectively. The value of ci for asymptotically flat met- 
rics is in this case 



Cl 



1024 \p, 



Pi 



-'1/2 



(84) 



and the relative difference to the Kerr quadrupole mo- 
ment is 



AQ _ 25pi 
<9 16pc 



(85) 



which diverges as 1/A. By scaling the input parameters 
Poc, Poc and luq in the numerical code used in the previous 
section as A, A^ and A^/^ respectively, these values are 
also obtained numerically in the limit A ^ 0. 



For Newtonian polytropes, p — Pc 



, where pc and 



Pc are the central density and pressure respectively, the 
integration of l|7^ gives the following equation for Ai 

^(fAi 4 dAi , , , , 1 

p = 2—^ + -—^ = k{air, - Ai)— (86) 
dr^ r dr 

with = Ai{ri) being the value of Ai at the zero 
pressure surface and 



Pc 7 



(87) 



The homogeneous version of (|77p then becomes 
d2/i2i/i . 2 (i/i2ih 



1 

^2 



r rfr 



(88) 



2 12 
^ (ai^i - Ai)^ + — ] h2ih = 0. 



From (El) it follows that 



d^A, 
dr^ I 



idAi 
dr 



-2/ri 



(89) 



on the zero pressure surface r — ri. Substitution of this 
and dSni) into ^ now gives 



Cl 



5 r^bj 



I,, {^/h2in - 2h 



48 M5 



(^^21. + 3/, 



(90) 



Let h2ih 7^ be the solution of ((88)) that makes the 
spacetime asymptotically flat. Substitution of /121 ~ 



oj^i.^r^ -f h2\h in (|4T|) then gives 



9i 



1 



9^21/1 + 3r 



C^^2U 



15r2 V'"'"' ' dr ""^Z^' 
from which the denominator of ([50)1 is 

'(i/i2i?i 



dr 



h2ih 



5^1/2 
3/I2I/1 



(91) 



(92) 



Note that /i2ift = would give a non-zero qi. 
With the identifications 7 = 1 + 1/n, 



-A<I> + airi , /l21/i 



A 



(93) 



in terms of the functions $ (that corresponds to the 
Newtonian potential for the non-rotating configuration) 
and 02 used in [7], and a rescaling of the r-coordinate 
r = \/2/k s, with k = fcA"^^, equations ([55]) and 
read as 



and 



ds2 



d2$ 

2 d(j)2 
s ds 



-? + $" = (94) 
s ds 



_ _ ) 02 = (95) 



respectively. Hence these are exactly the same as equa- 
tions (23) and (41) in the paper Q by Martin et.al. By 
rescaling the time coordinate to first order in A the cen- 
tral data for their function $ can be met and hence we 
may rely on their theorem on page 10. From their proof 
the following inequality holds for regular solutions 02 of 



0gsi) _1 

02(si) Si 



<0, 



(96) 



with 02 — d(j)2/ds and si 
dh2ih/dr 



'-21/1 



'fc/2ri, or in our notation 
2 



< 0. 



(97) 



That 02 (si) and hence also /i2ih(^i) 7^ is guaran- 
teed by their lemma. But (|57)) is actually exactly one of 
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the factors in the numerator of (|90|) and hence we con- 
clude that Ci 7^ 0. 

Substitution of (l92|) in ([90|) now gives that ci > iff 
h2ih\r=ri > Oj or equivalently if the fluid ball is oblate 
in shape. Hence AQ/Q = l&M^ci/ba? > for oblate 
Newtonian polytropes. 

VII. CONCLUSIONS 

It seems that generically the modulus of the 
quadrupole moment for slowly and rigidly rotating per- 
fect fluid balls is larger than that of the Kerr metric. 
Also, typically it becomes large for small central pres- 
sures. For the incompressible case the ratio AQ/Q is 
a monotonically decreasing function of central pressure. 
In the limit, when central pressure goes to infinity, it 
approaches a small, but positive value. Newtonian poly- 
tropes show a slightly more complicated behaviour, and 
AQ/Q possesses a minimum for 7 < 2. However, these 
minimal values are always larger than 0.2. The relativis- 
tic polytropes behave in a similar way, and the values 
of the minima are even higher. For linear equations of 
state something interesting happens. When the constant 
di is less than —1, the configuration becomes infinite in 
extent for a finite central pressure. In this limit the zero 
pressure surface approaches the event horizon and AQ/Q 
goes to zero. Also the central regions of these fiuids ap- 
proach the anti-de Sitter spacetime. All of these solution 
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